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Abstract 

Background: The imbalance between Th17 and Treg cells has been studied in various diseases including allergic 
asthma but their roles have not been fully understood in the development of the late phase asthmatic response. 

Objectives: To determine changes in Th17 and Treg cell numbers between isolated early responders (ERs) and dual 
responders (DRs) undergoing allergen inhalation challenge. To identify gene expression profiles associated with 
Th17 and Treg cells. 

Methods: 14 participants (8 ERs and 6 DRs) with mild allergic asthma underwent allergen inhalation challenge. 
Peripheral blood was collected prior to and 2 hours post allergen challenge. DNA methylation analysis was used to 
quantifiy the relative frequencies of Th17, Tregs, total B cells, and total T cells. Gene expression from whole blood 
was measured using microarrays. Technical replication of selected genes was performed using nanoString nCounter 
Elements. 

Results: The Th177Treg ratio significantly increased in DRs compared to ERs post allergen challenge compared to 
pre-challenge. Genes significantly correlated to Th17 and Treg cell counts were inversely correlated with each other. 
Genes significantly correlated with Th177Treg ratio included the cluster of genes of the leukocyte receptor complex 
located on chromosome 19q 13.4. 

Conclusions: Th17/Treg imbalance post-challenge may contribute to the development of the late phase inflammatory 
phenotype. 

Keywords: Allergen inhalation challenge, Asthma, Asthmatic response, DNA methylation, Epigenetic cell counting, 
Peripheral blood, Th17/Treg ratio, nCounter Elements 



Introduction 

The imbalance between a proinflammatory T helper 17 
(Thl7) and a regulatory T (Treg) cell phenotype may play 
a crucial role in allergic airway inflammation [1], Experi- 
mental models have shown that Thl7 cells typically pro- 
mote neutrophilic inflammation, and also play important 
roles in airway hyperresponsiveness in concert with Th2 
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cells [2], In peripheral blood, Thl7 cell counts have been 
shown to be higher in subjects with allergic asthma com- 
pared to healthy controls [3,4]. The percentage of Thl7 
cells and IL-17 levels in peripheral blood have been shown 
to be significantly elevated 24 hours after allergen chal- 
lenge in dual responders compared to early responders or 
healthy controls [5]. On the other hand, Treg cells main- 
tain immune homeostasis and regulate immune responses 
to allergens by preventing excessive inflammatory responses 
[6]. Treg cells were originally identified as CD4 + CD25 + 
T cells with a function to suppress immune responses 
[7]. In order to identify Treg cells, FOXP3 expression 
as a specific marker has been used, however, it is also 
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expressed in activated non-suppressor T cells [8,9]. Low 
levels of the IL-7 receptor (CD 127) in combination with 
high expression of CD4 and CD25 can be used to isolate 
highly purified suppressive Tregs [10]. Recently, DNA 
methylation analysis of the Treg specific demethylation 
region (TSDR) within the FOXP3 locus has been used to 
enumerate Treg cells, [11] which have been shown to 
significantly correlate with CD4 + CD25 + CD127 l0 and 
CD4 + CD25 + CD127 lo FOXP3 + cells [12]. 

We have previously demonstrated that peripheral blood 
is a useful biological material with which to study changes 
in the blood transcriptome, proteome and metabolome of 
individuals with mild atopic asthma undergoing allergen 
inhalation challenge [13-16]. In the present study, we have 
used qPCR based DNA methylation analysis to estimate 
the number of Thl7 cells, Treg cells, T cells and B cells 
in peripheral blood of mild atopic asthmatics undergo- 
ing allergen inhalation challenge. In the same individ- 
uals, we also analysed gene expression profiles in whole 
blood using microarrays to identify genes correlated 
with each cell type. We hypothesized that changes in 
specific immune cell counts in peripheral blood would 
be associated with the allergen-induced late phase asth- 
matic response. 

Methods 

Study participants and allergen inhalation challenge 

The Institutional Review Boards of the participating institu- 
tions, University of British Columbia, McMaster University 
and Universite Laval, approved this study. Fourteen individ- 
uals were recruited as part of the AllerGen NCE Clinical 
Investigator Collaborative (Canada) and provided written 
informed consent to undergo an allergen inhalation chal- 
lenge. All participants were non-smokers, free of other lung 
diseases, and not pregnant. Diagnosis of asthma was based 
on the Global Initiative for Asthma criteria. Participants 
were diagnosed with mild allergic asthma, and only used 
intermittent short-acting bronchodilators for treatment 
of their asthma. Participants had a baseline FEVx > 70% 
of predicted, and the PC 2 o> provocative concentration of 
methacholine required to produce a 20% decrease in 
FEVi, was <16 mg/mL. 

Skin prick tests were used to determine allergies to 
cat, and the dose of cat allergen extract for inhalation. 
Methacholine and allergen challenges were conducted as 
triad visits. On the first and third day, participants under- 
went methacholine inhalation tests for assessments of air- 
way hyperresponsiveness (AHR) as described previously 
[17,18]. The allergen-induced shift (post/pre in PC 2 o) was 
evaluated as the change in AHR. On the second day par- 
ticipants underwent allergen inhalation challenge with ex- 
tracts of cat pelt or hair in doubling doses until a drop in 
FEVi of at least of 20% was achieved, then FEVi was mea- 
sured at regular intervals up to 7 hours post-challenge as 



described previously [19]. All participants developed an 
early response which resolved within 1-3 hours after chal- 
lenge. Participants that demonstrated a maximum drop in 
FEVi of greater than 15% between 3 to 7 hours after aller- 
gen inhalation were classified as dual responders (DRs). 
Participants having an FEV! drop of 10% that was still fall- 
ing at the end of the 7 hour observation period were cate- 
gorized as DRs if they also demonstrated a drop in PC 2 o 
(post compared to pre methacholine challenge). Partici- 
pants who showed neither a drop in FEVi > 15% between 
3 to 7 hours after challenge nor a decreased PC 2 o were 
classified as isolated early responders (ERs). 

Blood collection and isolation of RNA and DNA 

Peripheral blood was obtained immediately before and 
2 hours post-challenge in PAXgene Blood RNA tubes 
(PreAnalytiX, Qiagen/BD, Valencia, CA, USA) for RNA 
and in K2 EDTA Vacutainer tubes (BD, Franklin Lakes, NJ, 
USA) for buffy coat and complete blood count (CBC) 
measurements. Cellular RNA was purified from 2.5 mL 
of whole blood in PAXgene tubes according to the 
manufacturers protocols using the RNeasy Mini Kit 
(Qiagen, Chatsworth, CA, USA). Total DNA was isolated 
from whole blood or buffy coat from EDTA tubes using 
QIAamp DNA Blood Mini Kit (Qiagen) according to the 
manufacturer s protocol. 

Epigenetic cell counting of lymphocyte subsets using 
DNA methylation analysis 

Cell counting of lymphocyte subsets was performed by 
Epiontis (Berlin, Germany) using quantitative real-time 
PCR (qPCR) based DNA methylation analysis [20,21]. 
Briefly, bisulphite conversion [22] of genomic DNA result- 
ing in either CpG-variants (if DNA is methylated) or 
TpG -variants (if DNA is unmethylated) was performed. 
Each qPCR assay is specific for either the demethylated 
FOXP3 TSDR (for Tregs) or the demethylated CD3D/G 
(for T cells) or the demethylated IL17A (for Thl7 cells) or 
the B cell specific demethylated gene region (for B cells) 
templates, since the demethylated version of these regions 
have been shown to be exclusively present in Treg, T cells, 
Thl7 cells and B cells respectively. The other qPCR assay 
is specific for a control region within the GAPDH gene, a 
target that is demethylated in all cells. The GAPDH PCR 
assay serves as a "load control" as it estimates the number 
of "total cells" in a given sample. The percentage of Treg 
cells, T cells, Thl7 cells and B cells in a sample is calcu- 
lated as: 

Percentage of a particular cell-type = [Copy Equivalents 
as determined with the PCR assay targeting the cell- 
specific DNA target region (e.g. TpG TSDR )]/[Copy 
Equivalents as determined with the GAPDH qPCR 
assay (TpG GAPDH )] x [100] x [2 a) ]. 
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In the equation above, the "Copy Equivalents" as deter- 
mined by the cell-specific PCR assay corresponds to "Treg 
cells", or "T cells" or "Thl7 cells" or "B cells" copies, re- 
spectively. The "Copy Equivalents" as determined with the 
GAPDH PCR assay corresponds to the "total cell" copies, 
respectively. A factor of "100" is used to translate the re- 
sult into percentage of cells. 

a) Only for Tregs a factor of "2" is applied in the equa- 
tion to correct for the fact that each cell has two copies 
of the (demethylated) GAPDH gene but each Treg has 
just one copy of the demethylated FOXP3 gene. As 
FOXP3 is X-chromosomally located, each Treg holds 
exactly one copy of the demethylated FOXP3 gene. 
Tregs from male subjects hold one X chromosome on 
which the FOXP3 gene is demethylated. In contrast, 
each Treg from a female subject has two X chromo- 
somes (and thus two copies of the FOXP3 gene) but one 
X chromosome is inactivated (i.e. fully methylated) and 
it exists as a Barr body in the cell. 

Microarray gene expression assay 

Genome-wide expression profiling, labelling and array 
hybridization were performed using Affymetrix Human 
Gene 1.0 ST arrays (Affymetrix, Santa Clara, CA, USA). 
All microarray data has previously [16] been deposited 
into the Gene Expression Omnibus (GSE40240). All 
'CEL' files were normalized using the Robust Multiarray 
Average (RMA). 

nCounter Elements 

Technical replication of selected genes was performed 
using a new digital technology, nCounter Elements 
(NanoString, Seattle, USA). nCounter Elements allows 
users to combine nCounter Elements General Purpose 
Reagents (GPRs) with unlabelled probes that target spe- 
cific genes of interest (www.nanostring.com/elements/). 
100 ng of each RNA sample is added to the TagSet in 
hybridization buffer and incubated at 65°C for 16 hours. 
The TagSet consists of a reporter tag and capture tag that 
hybridize to the user designed gene-specific probe A and 
probe B complex. Automated processing per cartridge 
on the PrepStation (high sensitivity protocol) occurs for 
3 hours. After a 2.5 hour scan per cartridge, counts are 
acquired from the GEN2 Digital Analyzer. Details re- 
garding data normalization can be found in the supple- 
mentary material. 

Statistical and bioinformatics analysis 

Linear models were used to test the association between 
immune cell frequencies and cell-specific gene expression 
profiles. Cell counts and all combinations of cell ratios (T, B, 
Treg and Thl7) were compared using linear regression 
models. All microarray data were analysed using the lin- 
ear models for microarrays (limma) R-library [23]. The 



Benjamini-Hochberg false discovery rate (FDR) was 
used to correct for multiple testing. Partial least squares 
(PLS), from the mixOmics R-library [24] was used to 
identify the relationship between cell-specific gene lists. 
Statistical analyses were performed in the statistical 
computing program R version 3.0.1 [25]. 

To test for the enrichment of gene lists, GeneGo network 
analysis was performed using MetaCore from Thomson 
Reuters. Network analyses were performed on gene lists 
created by ranking genes by the scores which rank the 
subnetworks to saturation with the objects from the initial 
gene list. 

Results 

Participant characteristics 

The 14 participants were classified into eight isolated 
early responders (ERs) and six dual responders (DRs), as 
shown in Table 1. The mean drop in FEVi during the late 
phase in DRs (21.3 ± 3.2) was 4 times greater (p < 0.05) 
compared to ERs (5.1 ± 1.4). Table 1 also shows that all 
participants exhibited an immediate drop in FEV! of 
greater than 20%. 

Correlation between immune cell frequency and 
cell-specific gene expression 

Sum of the T cell and B cell frequencies obtained using 
the methylation assays strongly correlate (Spearman r = 0.95) 
with the lymphocyte frequency obtained using a hematolyzer 
(Additional file 1: Figure SI). T cell, B cell and Thl7 cell 
counts were significantly positively correlated with the 
genes targeted in epigenetic cell counting in both the 
microarray (Figure 1; top row) and nanoString (Figure 1; 
bottom row) platforms. Treg cell counts were not corre- 
lated with FOXP3 gene expression measured using micro- 
arrays but was significantly correlated using nanoString, 
suggesting greater sensitivity of the platform (Figure 1, red 
points). 

Th17 to Treg ratio discriminates early from dual 
responders after challenge 

Allergen inhalation did not significantly change T cell, B 
cell, Treg cell and Thl7 cell counts in either ERs or DRs. 
In addition, comparing the change in cell counts in ERs 
with the change in cell counts in DRs (AER vs. ADR) no 
significant cell- types were identified (Table 2A). Next, the 
ratios between different cell-types were analyzed (Table 2B). 
Table 2B shows that the Thl7/Treg ratio significantly 
(p = 0.03) increased in DRs compared to ERs, from pre 
to post challenge. Figure 2 shows that the Thl7/Treg 
ratio did not change from pre to post challenge in ERs 
(net change = 0.006 ± 0.09), whereas the Thl7/Treg ratio 
increased in DRs (net change = 0.28 ± 0.03). 
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Table 1 Participant demographics 



Participant 
ID 


Age (year) 


Sex (M:F) 


Allergen 


Pre [PC 20 ] 
(mg/mL) 


Post [PC 20 ] 
(mg/mL) 


Allergen-induced 
shift b 


% fall in FEV1 
Early Late 


ER 


















1 


28 


F 


Cat Pelt 


12.8 


ND 


ND 


20.3 


4.8 


2 


34 


F 


Cat Pelt 


2.8 


6.1 


2.3 


21 


1.5 


3 


27 


M 


Cat Pelt 


4.5 


1.8 


0.39 


34.4 


0 


4 


42 


F 


Cat Hair 


5.3 


8.6 


1.6 


42.1 


11.1 


5 


29 


F 


Cat Pelt 


0.4 


ND 


ND 


44.3 


0 


6 


31 


M 


Cat Pelt 


11.8 


16 


1.4 


24.2 


7.5 


7 


28 


F 


Cat Hair 


9.4 


16 


1.7 


27.1 


7.1 


8 


42 


M 


Cat 


0.1 


ND 


ND 


23 


9 


Mean±SE 


32.6 ± 2.2 


3:5 




2.8 a 


7.5 a 


1.5 ±0.3 


29.6 ± 3.2 


5.1 ±1.4 


DR 


















9 


23 


F 


Cat Hair 


0.3 


0.2 


0.60 


38.9 


31.8 


10 


26 


F 


Cat Hair 


5.1 


1.5 


0.30 


31.4 


14.9 


11 


49 


F 


Cat Hair 


3.6 


1.0 


0.27 


25.3 


12.6 


12 


26 


M 


Cat Hair 


0.9 


1.0 


1.1 


31.5 


15.6 


13 


27 


F 


Cat Pelt 


0.6 


0.1 


0.19 


48.3 


25.8 


14 


52 


F 


Cat 


ND 


ND 


ND 


33 


27 


Mean±SE 


33.8 ±3.3 


1:5 




1.3 a 


0.5 a 


0.5 ±0.1 


34.7 ± 3.2 


21.3±3.2 C 



ND - Not determined; all participants were challenged with cat allergen. 
a geometric mean (PC 20 values are measured on a log scale). 

[PC2o]post/[PC 2 o] pre- 

c p < 0.05 versus ER group. 



CD3D vs T-cell 



CD3G vs T-cell 



CD79A vs B-cell 



CD79B vs B-cell 



FOXP3 vs Treg 



IL17A vs Th17 




Cell frequency 

Figure 1 Scatter plots of immune cells quantified using DNA methylation analysis with the corresponding cell-specific gene expressions 
profiles, x-axis: relative cell-type frequencies of T, B, Treg and Thl 7 cells in whole blood; y-axis: a) top row: gene {CD3D, CD3G, CD79A, CD79A, 
FOXP3, and IL17A) expression intensities measured using microarrays and b) bottom row: gene expression counts measured using nCounter 
Elements from nanoString. 
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Table 2 Comparing immune cell-frequencies and cell/cell ratios between early and dual responders after allergen 
challenge 



A. Comparing immune cell-frequencies between ERs and DRs after allergen challenge 




Fold-change (post-pre) in early responders 


Fold-change (post-pre) in dual responders 


AER vs. ADR 




Mean ± SEM 




Mean ± SEM 


p-value 


Tcell 


-2.20 ± 1 .40 




-4.26 ± 2.39 


0.45 


Bcell 


-0.20 ±0.1 7 




-0.11 ±0.26 


0.77 


Th17 cell 


-0.22 ±0.21 




0.14 + 0.16 


0.22 


Treg cell 


-0.1 2 ±0.06 




-0.42 ±0.1 7 


0.10 


B. Comparing immune cell/cell ratios between ERs and DRs after allergen challenge 




Fold-change (post-pre) in early responders 


Fold-change (post-pre) in dual responders 


AER vs. ADR 




Mean ± SEM 




Mean ± SEM 


p-value 


T/B 


-0.59 ± 0.70 




-0.74 ± 0.43 


0.87 


Th17/T 


0.002 ± 0.006 




0.009 ± 0.002 


0.31 


Treg^ 


-0.0008 ± 0.005 




-0.006 ± 0.002 


0.38 


BATh17 


-0.009 ±0.10 




-0.37 ±0.21 


0.13 


B^reg 


-0.002 ± 0.22 




0.42 ±0.1 8 


0.18 


Th17^reg 


0.006 ± 0.09 




0.28 ± 0.03 


0.03 



Genes associated with Th17 and Treg cells 

A multiple linear regression model (limma) was used to 
identify genes whose expression levels correlated with 
the frequencies of specific cell-types independent of 
changes in the frequencies of other cell-types (gene ex- 
pression ~ Thl7 + Treg + B-cells + other T-cells, where 
other T-cells = overall T-cells minus Thl7 and Treg). 
10 (99) genes were positively correlated with Thl7 (Treg) 
cells at an FDR of 10%, with no overlapping genes 



between the two lists. Thl7 genes included KIR2DS2, 
TAGLN, C14orf37, KRTAP13-3, SAP30, KIR2DS4, LAIR2, 
FLJ30679, RORC and KIR2DL2. The 99 Treg genes were 
enriched (FDR = 5%) for 27 pathways including many 
relevant regulatory pathways such as IL-2 regulation of 
translation, Regulation of telomere length and cellular 
immortalization, Regulation of T cell function by CTLA-4 
(Additional file 1: Figure S2). Partial Least Squares (PLS) 
was used to determine the correlation between the set of 



2.0 



K 1.6 



Early Responders 




1.5 



1.0 



Dual Responders 




post 



Time 



pre 



post 



Time 



Change in Th17/Treg ratio 



- - -e- 



n 




Net change in ERs 
Mean=0.006, SEM=0.09 

Net change in DRs 
Mean=0.28, SEM=0.03 

ER vs. DR; p=0.03 



ER DR 

Response 



Figure 2 Change in the Th17/Treg ratio in early and dual responders from pre to post challenge. Th17^reg ratio in ERs (left panel) and 
DRs (middle panel) at pre and post-challenge. The change in the Th17^reg ratio (post-pre) in ERs and DRs (right panel). Solid black points 
depict the median value of the data for each boxplot. 
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10 Thl7 genes and the set of 99 Treg genes. Figure 3 de- 
picts the results of PLS using a correlation circle (see 
Gonzalez et al [26] for complete details on graphical 
outputs of PLS). Vectors drawn from the origin to each 
of the points (genes) allows one to determine the relation- 
ship between genes: 1) if the angle between two vectors is 
less than 90°, there exists a positive correlation between 
the two genes, 2) if the angle between two vectors is 
greater than 90°, there exists a negative correlation be- 
tween the two genes, and 3) if the angle between two vec- 
tors is equal to 90°, the correlation between the two genes 
is zero. Figure 3 shows that the Thl7 genes were inversely 
correlated with Treg genes (angle greater than 90°). 

Genes significantly correlated with the Th17/Treg ratio 

To investigate the relationship of Thl7/Treg ratio and 
gene expression profiles, we identified correlated genes 
in the entire sample set. We identified 13 genes significantly 
correlated with Thl7/Treg ratio using limma (FDR = 5%, 
Table 3). Interestingly, 7 genes (KIR3DL1, LAIR2, KIR2DS2, 
KIR2DL2, CD226, KIR2DS4, KIR2DS1) belong to the 
leukocyte receptor complex (LRC) located on chromo- 
some 19ql3.4, and were shown to be positively corre- 
lated except CD226. However, of the four genes profiled 

c s \ 











/,FGFR1QP 


ThtX 


// ZNF33A \ 








2S6Bjb8£ >P7 / \ 


\ / TAGLN y 






j KRTAP1 3-3 




- T^SjjZ&y2 Ar 1 Lr 


I ! LAIR2 




^WWPJpSbskiaaixj? / 


J \ RORC / 
/ \ C14orf37 
/ \ SAP30 
/ \U30679 I 












1 1 


I I 



-1.0 -0.5 0.0 0.5 1.0 



Comp 1 

Figure 3 Correlation circle depicting the strength of correlation 
between Treg genes (red) and Th17 genes (blue) with their 
respective latent variables (Comp 1 and Comp 2). The Treg 
genes (red) show a strong negative correlation with the Thl 7 genes 
(blue). Vectors drawn from the origin to each of the points (genes) 
allows one to determine the relationship between genes: 1) if the 
angle between two vectors is less than 90°, there exists a positive 
correlation between the two genes, 2) if the angle between two 
vectors is greater than 90°, there exists a negative correlation 
between the two genes, and 3) if the angle between two vectors is 
equal to 90°, the correlation between the two genes is zero. 

V ) 



Table 3 Genes significantly correlated to Th17/Treg ratio 
in Pearson tests (FDR <0.05) 



Gene 


X 


p value 


FDR 


TAGLN 


0.78 


7.59E-07 


0.02 


C14orf37 


0.77 


3.16E-06 


0.02 


KIR3DL 1 * 


0.75 


3.82E-06 


0.02 


LAIR2* 


0.75 


4.80E-06 


0.02 


CDCP1 


0.74 


5.77E-06 


0.02 


KIR2DS2* 


0.74 


6.19E-06 


0.02 


SAP30 


0.74 


7.09E-06 


0.02 


KIR2DL2* 


0.73 


8.51 E-06 


0.02 


CD226* 


-0.73 


9.83E-06 


0.02 


ZNF286B 


-0.73 


9.27E-06 


0.02 


KRTAP13-3 


0.72 


1.57E-05 


0.03 


KIR2DS4* 


0.71 


2.38E-05 


0.043 


KIR2DS 1 * 


0.70 


2.86E-05 


0.048 



*genes belonging to leukocyte receptor complex (LRC). 



using nanoString, only CD226 and KIR2DS4 success- 
fully replicated (Figure 4). The top-listed transcriptional 
network in GeneGo network analysis for the 13 significant 
genes included regulatory functions in immune responses 
(Additional file 1: Table SI). 

Significantly different genes between ERs and DRs in 
Thl 7 or Treg 

In a secondary analysis, we also analysed gene-cell correla- 
tions that significantly differed between early and dual re- 
sponded, irrespective of allergen exposure, using limma. 
In GeneGo network analysis, genes differentially associ- 
ated with Thl7 (165 genes differentially associated with 
Thl7 between ERs and DRs, p < 0.01) were enriched for 
immunological processes including immunoglobulin me- 
diated immune response and adaptive immune response. 
Genes differentially associated with Treg between ER and 
DR (554 genes, p < 0.01) were enriched for immune pro- 
cesses. Although the genes differentially associated with 
Thl7 cells between ERs and DRs did not achieve a strin- 
gent threshold of FDR, the top three genes, SIOOB, MILR1 
and CHI3L1 (p-value < 0.001, FDR = 0.79, Additional file 
1: Figure S3), have previously been reported to be involved 
in allergy or asthma [27-29]. Additional file 1: Figure S3 
shows that all three genes were differentially correlated 
with Thl7 cell counts with respect to the response class 
using both the microarray and nanoString platforms. 

Discussion 

Although Thl7 and Treg cells arise from a common 
precursor cell [30] they have opposing inflammatory 
roles which has been demonstrated in the context of 
autoimmune disease [31], infection [32], and recently 
allergic airway inflammation [1]. In the present study 
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CD226 



R = 0.52 
Estimate = -0.48 
,SEM = 0.088 



R =0.39 
Estimate = -0.5 
SEM = 0.12 
p<0.001 
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2.0 



0.5 1.0 
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2.0 
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• 
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• 


SEM = 0.081 
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• • 
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• 


' • • • 
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R 2 = 0.43 
Estimate = 1.7 
SEM = 0.36 
p<0.001 



KIR3DL1 



fT = 0.55 
Estimate = 0.59 
SEM = 0.1 



R =-0.013 
Estimate = -0.22 
SEM = 0.28 i 
p=0.09 



0.5 

Th17/Treg ratio 



1.0 

i_ 



LAIR2 



R =0.54 
Estimate = 0.44 
SEM = 0.077 



K = 0.075 
Estimate = -0.23 
SEM = 0.13 # 
p=0.43 



"T" 



Figure 4 Scatter plots of genes significantly correlated with the Th17/Treg ratio, x-axis: Th17/Treg ratio; y-axis: a) top row: gene {CD226, 
KIR2DS4, KIR3DL1, and LAIR2) expression intensities measured using microarrays and b) bottom row: gene expression counts measured using 
nCounter Elements from nanoString. 



we demonstrate a potential Thl7/Treg homeostatic im- 
balance using peripheral blood of isolated early and dual 
asthmatic responders (ERs and DRs) undergoing allergen 
inhalation challenge. 

DNA methylation analysis used to enumerate various 
immune cells revealed good correlation with the cell- 
specific gene expression profiles as measured using micro- 
arrays. Technical replication using nCounter Elements 
from nanoString, a more sensitive platform indicated that 
FOXP3 expression was indeed correlated with Treg cell 
counts. As a marker for human Tregs, however, FOXP3 
expression is of doubtful value, due to its transient expres- 
sion in activated non-regulatory effector T cells [21]. In 
addition, other cell-surface markers such as CD 127 or 
CD45RA have been used to isolate FOXP3 + Treg cell pop- 
ulations with high efficiency [33,34]. Epigenetic enumer- 
ation of Treg cells in the present study has been shown to 
positively correlate with CD4 + CD25 + CD127 lG , and CD4 + 
CD25 + CD127 lo FOXP3 + [12] and thus are truly repre- 
sentative of suppressive Tregs. 

The percentage of Treg cells did not significantly change 
in either ERs (-0.12 ± 0.06; p = 0.11) or DRs (-0.42 ± 0.17; 
p = 0.054), two hours post-challenge. Previous studies have 
also not shown significant changes in Treg cells in periph- 
eral blood in DRs undergoing allergen inhalation chal- 
lenge [35,36]. This may be due to many factors such as 
the time of the post-challenge blood draw, the cell-surface 
markers used to isolate the Treg cells as well as the small 
sample sizes (n = 6-11) used in these studies. Similarly, the 
percentage of Thl7 cells also did not significantly change 
in ERs (-0.22 ± 0.21; p = 0.30) or DRs (0.14 ± 0.16; p = 0.44), 



after allergen challenge. Thl7 cells have been shown to be 
increased 7 and 24 hours post-challenge in both ERs and 
DRs and the increase in DRs was greater than in ERs 
24 hours post challenge [5]. Thl7 cells as well as the con- 
centrations of IL-17 and IL-22, have also been shown to 
be increased with the severity of allergic asthma [37]. 
Genes significantly positively correlated with Thl7 cells 
included RORC, the transcription factor involved in Thl7 
differentiation, whereas genes significantly positively cor- 
related with Treg genes was enriched for regulatory func- 
tions. Furthermore, Thl7 and Treg cell associated genes 
were inversely correlated with each other, further implicat- 
ing the phenotypic roles of these cell-types in allergic 
asthma. 

Although neither cell-type significantly changed pre to 
post challenge, the change in the Thl7/Treg ratio from 
pre to post challenge significantly (p = 0.03) differed be- 
tween ERs and DRs. The Thl7/Treg ratio increased in 
DRs whereas little change occurs in ERs after challenge. 
The increase in the Thl7/Treg ratio in DRs is driven by 
an increase in the number of Thl7 cells (0.14 ±0.16) 
and a decrease in the number of Treg cells (-0.42 ± 
0.17) due to allergen exposure. A possible mechanism of 
Thl7/Treg imbalance was suggested by the genes that 
were correlated with Thl7/Treg ratio. LRC on chromo- 
some 19ql3.4 encodes immunoglobulin super family re- 
ceptors including killer immunoglobulin like receptors 
(KIRs) expressed on hematopoietic cells. Almost all LRC 
significant genes were positively correlated to Thl7/Treg, 
whilst CD226 is the only LRC gene negatively correlated. 
A previous study on differential expression of LRC genes 
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revealed that KIRs and inhibitory receptor ILT2/LIR1 
were expressed in activated T cells and that KIR levels in 
T cells are associated with resistance to activation-induced 
cell death [38]. These may suggest a new hypothesis that 
LRC gene expression patterns might be related to Thl7/ 
Treg ratio and involved in immune responses to inhaled 
allergen in asthmatics. 

The statistical interaction analyses suggested differences 
in gene expression profiles in Thl7 or Tregs between ERs 
and DRs. Interestingly, top-listed differentially expressed 
Thl7-associated genes SIOOB, MILR1 and CHI3L1 have 
been reported to play roles in allergy and asthma. S100B + 
lymphocytes in blood have been reported to consist of 
two subtypes; a cytotoxic T cell and a NK subtype [27] . In 
connection with the significant correlations between Thl7 
and KIR family, Thl7 measured by epigenetic cell count- 
ing for IL1 7 A might be related to other types of immune 
cells. This is supported by reports showing that IL-17 
genes are expressed in non-CD4 + T cells such as yd T 
cells, NK cells and Type 3 innate lymphoid cells, suggest- 
ing that innate immunity might be responsible for initiat- 
ing this type of inflammation commonly associated with 
Thl7 immunity [39,40]. Further studies are needed to 
clarify the disparity between true Thl7 and IL17A- 
demethylated cells. MILR1 is the gene for allergin-1 pro- 
tein, which was recently identified to play an inhibitory 
role in mast cell functions [28]. Polymorphisms in CHI3L1 
as well as the concentration of its corresponding protein 
YKL-40 in serum has been associated with asthma and 
pulmonary function [29]. Our findings suggest that Thl7 
cell gene expression profiles are divergent between asth- 
matic responses and that these profiles might be related to 
immune mechanisms. 

A limitation of this study is its small sample size, 
which reduces the statistical power in identifying true 
positives. Therefore we deemed a technical validation 
using a highly sensitive platform appropriate for this 
study. Independent replication will be important as part 
of future studies with larger sample sizes. Another limi- 
tation of the present study is that only a limited number 
of cell-types were studied using DNA methylation ana- 
lysis, whereas quantification of a wide array of cell-types 
such as Thl, Th2, and Th9 cells would provide deeper 
biological insights into the mechanisms of allergic asth- 
matic responses. DNA methylation based qPCR assays for 
these cell-types will allow for tissue samples to remain un- 
perturbed, and additional sources of variability, such as 
those observed in fluorescence activated cell sorting to be 
avoided. 

The careful phenotyping of our participants, together 
with innovative epigenetic- and gene expression-based 
methodologies, have nevertheless revealed interesting di- 
rections for further investigations using large sample 
sizes and different allergens. 
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